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Abstract 

The interaction energy of the two-dimensional electron system in the region of fractional 
quantum Hall effect is considered within the Chern-Simons composite fermion approach. 
In the limit when Coulomb interaction is very small comparing to the cyclotron energy 
the RPA results are obtained for the fillings v = 1/3, 1/5, 2/3, 2/5, 3/7 and compared 
with the exact diagonalization results for small systems (extrapolated for infinite sys- 
tems). They show very poor agreement suggesting the need for looking for alternative 
approaches. 

PACS: 73.40.H, 71.10.P 

keywords: fractional quantum Hall effect, composite fermions 



1 



1. Introduction 



One of the aims of tlieoretical studies of fractional-quantum-Hall-effect (FQHE) sys- 
tems [|l], ^ is to determine the value of the interaction energy of the system at various 
fractional fillings of the lowest Landau level. The exact numerical results can be found 
for few particle systems ^, Q and they agree very well with the predictions of the 
trial wave function approach of Laughlin and Jain |I|, ^, |^ (Laughlin and Jain wave 
functions are proposed for the case when Coulomb interaction is very small comparing 
to the cyclotron energy and higher Landau levels can be omitted, it is also the case in 
related numerical work). The many-body theory for such systems is formulated within 
the Chern-Simons gauge theory Q which introduces the gauge field mapping fermions 
into fermions (so-called composite fermions obtained by attaching an even number of flux 
quanta to each electron). There is no small parameter in the Chern-Simons composite 
fermion theory |10|, |Tl], |l^ (in contrast to the case of anyons in the fermion limit |1I3| ). 
Nevertheless, the Chern-Simons theory gives good predictions of the transport proper- 
ties of the system TH]. However, in contrast to the trial wave function approach, the 
attempts to get the ground state energy of the system within the Chern-Simons theory 



are not very successful |]T2|, |T^. In calculations of energy gaps modified approach have 
been used in which the relation between finite size calculations and the Chern-Simons 



results is assumed M, 11, O. In this paper we calculate the interaction energy of the 



system in the Chern-Simons approach within the RPA for several values of the filling 
fraction (|, |, |, |, |) and compare them with the exact diagonalization results for few 
electron systems (extrapolated for infinite systems). 

The composite fermion (CF) transformation consists in attaching an even number 
(2p) of flux quanta to each electron. Such point fluxes do not change the statistics 
of composite particles, however, they allow to treat the 2D system of electrons in a 
strong magnetic field in a close analogy to the treatment in a weak magnetic field. It 
is motivated by the mean field approach when the sum of point fluxes is replaced by an 
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average flux. The corresponding average field is B'~^^ ^ = —2p^p (fluxes opposite to 
the external flux, p is an integer, p - density) and the effective field acting on electrons 
is reduced to B'f^ = fi- + 5^^-^. In a close analogy to the quantum Hall effect one 
predicts similar effect when n Landau levels are completely filled in the effective field, i.e. 
B^ff = Hence, 5^^' = 2p!l±i^p ^hich means that the "real" lowest Landau level 

is filled in the fraction v = 2pn+i ~ 2pn-i ^^en the effective field is opposite to the 
external one). It is interesting to notice that the Laughlin fractions (of the form 1/m, 
m - odd) can be represented in two different ways. We can add that the Chern-Simons 
theory results should be independent of the actual Chern-Simons parameter 2p (if even) 
p!6| . Hence, we can use the value of p which is most suitable in a given problem. In 



practice, we use the value of p which gives fully filled Landau levels in the effective field 
(treatment of such systems is well known). Nevertheless, whatever value oi p is taken 
the results should be the same. 

The Hamiltonian of the two-dimensional system of electrons in an external magnetic 
field 

H= /t;2r^+(r) — (p + -A""(r))2^(r) 
J 2m c 

+ - f f dhdh'-^+(r)^+(r')—^—--^(rW(r') (1) 
2 J J e|r — r'l 

(e - dielectric constant) can be rewritten in the following way: 

H = [ dh^+{r) — (p + -A^"(r) + -A^^-^(r))2^(r) + 
J 2m c c 

where 

A^^'-'ir) = -2p^ [ dhej^^pir) , (3) 
e J I r — r I ^ 

p(r) = \I'+(r)\E'(r). The Hamiltonian H can be separated into two parts: H = Hq + Hint 
where 

Ho = J dh^^{r)^ip + -^A^f{r)f^{r) (4) 
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is treated as the unperturbed term {B^-f = V x A''-^ = V x ( A''^+ A'^''"'^) = B''^+B'^^~^ , 
QCh~s jg fQ^j^fi \jy averaging point fluxes - putting the average density p in (3)). Hint is 



the interaction Hamiltonian [[13], |Tj] : 

H.nt = \j j rf2rd^r'*+(r)vl/+(r')^^:^*(r)vl/(r') + H, + , (5) 

where 

H, = -2p- f f dhdh'^+{v){p^ + ^-A:f)^{v)eJ^^^{p{r') - p), (6) 



m 



In this paper we consider the case when B'^^^ = ^^p, i.e. in the unperturbed state 
one has n completely filled Landau levels (the effective filling u* = n). The first step 
in calculating the ground state interaction energy is the Hartree-Fock approximation. 
Considering the Coulomb interaction one finds the Hartree-Fock (H-F) contribution to 
be (we assume the presence of the positive background): 



E 



H-F 



/ [E^°(^0]'exp(--Orfr (8) 



^0 ' 



where a^-^-^ = y/j^irr is the effective magnetic length (cg'^'^ = \f^^al^ = ^J2pn ± la; 
Oq^ = \J ^^), -Z^r - Laguerre polynomials, - number of particles. The H-F results 
are presented in Table I for several filling fractions and compared with "exact" results 
(exact diagonalization results extrapolated for infinite systems) . The difference between 
"exact" and Hartree-Fock results (correlation energy) increases with the decrease of the 
fraction, for the 1/3 state is of order of 10% (of the exact value). It seems that a higher 
order approximation, eg. the RPA, will give a better agreement. In the following we 
consider the correlation energy within the RPA, assuming that the separation between 
Landau levels is much larger than Coulomb interaction between particles (as it is the 
case in exact diagonalization methods we refer to). 
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2. Correlation energy 



The correlation energy can be defined as follows: 

"1 dX 



/ — (< XHint >A - < Aifi^t >o) . (9) 
Jo A 



The expression for the correlation energy in the RPA (three-body contributions are 



omitted) has the form p8| , [12 1 

Er^ = -hL^J ^^°°^/^'^Imtr(Ang))[/^f^^(q,o.)-/^o(q,^)] (10) 

where D^^"^ is the correlation function of effective field currents (L^ is the area of the 
system) : 

^2f^(rt,rY) = ~< T[firt),fir't')] > (11) 
given within the random-phase approximation (with the coupling constant A): 

Dr^(q,a;) = [I-XDoiq,u;)V{q)]-'Do{q,u;). (12) 

The current densities are defined as: 

jW = ^E{P. + ^Af,5(r-rj)} (13) 

where braces denote an anticommutator, j is the vector part of with fi = 0,x, y. We 
define j° as density fluctuations: j° = J2j ^(r — Tj) — p. The interaction matrix V is 
obtained from the Hamiltonian Hint (dropping out three-body terms). We choose q = gx 
and the Coulomb gauge which reduces the problem to 2 x 2 Dj^^^ matrix {fi = 0,y) p. 
Taking ul^^ = and a^-^-^ to be frequency and length units, respectively one finds 

(^ = 1): 



^(q) 



^ v{q) \ Aprc I 2pn -iq 



V 



0/ 



(14) 



where v{q) is the Fourier transform of the Coulomb potential {v{q) = in standard 
units). 



5 



Let us assume the correspondence between the two energy scales, one related to the 
separation between Landau levels {hu'^^), one to the strength of the Coulomb interaction 

2 

(^fsr)- We introduce the dimensionless parameter: 

1 

(15) 



* ~ e// f- eff 
eOg huJc 

which shows the strength of the Coulomb interaction with respect to the separation 
between effective Landau levels (r^^ = ■ = ^2pn±i ^^'^' '^^^ considers the 
system of electrons, the limit ■ — > corresponds to the case when particle-hole 
excitations (electron excited into a higher Landau level) are negligible (hence, in exact 
diagonalization studies higher Landau levels can be omitted). When applying the Chern- 
Simons picture, however, gauge interactions are always of order of hu^^, and particle-hole 
excitations (CF excited into an empty Landau level and a CF hole in a filled level) have 
to be considered. We have 



1 \ iq 



The correlation function Dq is [|T9| : 



Tl 

^o(q,^) = ^ 




(17) 



where 



-X oo n-1 m — ] /' 
V _ \^ \^ ^ ^_ m-l-l^jm-U 

' n hnhW-irn-l-^ri)-rr^'' ^""^ ^^^^ 

X [(m - / - x)L["-'(x) + 2x F (18) 



and X = %r. Then one obtains: 



27cdet 



where det = det{I - D^V) = (1 - 2pnSi)2 - (2pn)2So(l + S2) - nr.gSo, 

= El - 2pnT,l + 2pnSoS2, Sp = {2pnfT,l + S2 - (2pn)2SoS2 + gnr,(S2 - S0S2). 



(19) 



Collective modes are determined by the poles of the correlation function D^^^. In 
Figures 1-2 we plot the results for v = 1 (the direct result and the p = 1 CF approach 
result) and similar results for u = 1/3 are presented in Figures 3-4, in Figure 5 the 3/7 
case is presented. 

The RPA correlation energy will be found using the dispersion relation of collective 
modes. In units of — ^ the correlation energy can be expressed as follows [|l8l {-ttt = 
TshoJl^^ ): 

^RPA^ / qdq — Im(lnrfet + tr(V(g)Do)) (20) 

2nrs Jo Jo tt 

which equals 

J\J POO pca fj^^j 

^RPA ^ / / — iui(\ndet + 2pn(2pnSo + 2Si) + nr.gSo). (21) 

ZriTs Jo Jo TT 

We have: 

dx / — ImSo(x, uj) = --Y.- + -{Sn - 1) (22) 

{Sn = Ej=i j)- This term is divergent but combined with other terms in the integral 
(21) has to give a finite value. Additionally 

/ dx —\mLi{x,uo) = {]. (23) 
JO JO vr 



The last integral in (21) 



^r,^,,r;^toE, (24) 

I Jo Jo Tl 



will be calculated separately for different v* = n. 



3. Results 



To calculate the correlation energy (21) one needs to know the zeros of the determi- 
nant det (collective modes). For n = 1 (the effective filling z/* = 1 - Laughlin fractions) 
we have an infinite set of modes with shortwavelength behaviour like uJmiQ oo) = m 
- Figs. 1-4 |]TD[. It can be shown that |TB[: 



f°° duj °° °° 

/ — Im In det = ^ {cUm " '"^) = X! (25) 

•^0 ^ m=l m=l 



The integral ( p^ equals: 

duj roo ^ e-^x^-i 



/ q^dq / — ImSo = / E , = / v^(l - e"^) • (26) 

7o JO TT Jo m! 7o 

Again, this integral is divergent, but combined with (25) and (22) will give a finite result. 
We write: 

q-dq ^ImSo = 4V2^E^n(2^-l) (27) 



/ — Im In det = Acui + V ( Acu„ + Acj+ ) (29) 



and then (in units of -^fj) 

= — y ( / ^i^m{x)dx - {2pf— - -r.v^ft ^^) . (28) 
2r, ^ ^ ^ 2m 2 2z ^ ^ ^ 

* m=l 1=1 

In the case when n > 2 every root (of collective modes) higher than the first is splitted 
into two |jl8| (for m > 1, uj^{q oo) = m = uj^{q oo) - Figure 5). One has 

'■"^ dw^ 

Im In det = Aui + 

m=2 

2 

and the correlation energy for n = 2 is given by (in units of — 177) 
-IT = V A.;i(x)dx - 8p') + ^ E [/ + ^^M)dx - [Apr—] + ^ 

^ fT^m A ^ ^ (2m + 3)(2m + l) ^ ™ 2^ - 1 

We have also found the expression for E^^^ for n = 3 (applied for = 3/7). 

The value of interaction energy related to Coulomb interaction (in the limit of 
Tg — > 0) is: 

= .H^^V N ) ' 

A^C stands for the case with no Coulomb interaction WM. . For n = 1 one has 



A rpRPA 1 00 „ 1 ™ 9,' _ 1 

^ ^ - ^-''(^))^^ - ^^^^ n ^) • (32) 



■» m=l j=l 

The main problem in calculating (32) is the calculation of the integrals and the conver- 
gence in summation over m. The integrals in (32) (and similar for n = 2 and n = 3) 
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have been calculated numerically using /c-point Gauss-Laquerre integration. It was veri- 
fied that the summation over m converges well and the sums have been truncated at 2k 
terms. In order to find the limit — > we considered small values of r^. It appears 
that for r, of order of 10"^ - 10"^ the RPA energies become practically independent of 
and for that range the results are given in Tables II-IV. 

2 

In Table II we present the RPA results in respective units of — 177. It can be seen that 
the results for the same value of p and are very close one to the other. They are not 
the same for the same filling (as they should be). We observe rather strong dependence 
on p (at a given effective filling u*). 

In Table III and Table IV the RPA results are compared with the exact diagonaliza- 
tion results. The best agreement with the "exact" values is found for the series 3/7, 2/5, 
1/3 (fractions going down from 1/2). The 1/3 state gives the best result, the difference 
between the RPA and the exact values is of order of 20% (of the exact result). 

An interesting example is the system of electrons at the real filling u = 1. In the limit 
r^^ — y the RPA correlation energy is zero. The qualitative difference is found within 
the RPA for the CF approach. Performing the CF transformation one finds the system 
at the effective filling u* = —1 (the effective field is opposite to the external one). We 
plot collective modes for the two descriptions in Figures 1-2. Calculating the RPA result 
(as it is described above) we find a finite CF result for the correlation energy (r^ — > 0). 

Similar situation can be found at the filling u = 1/3. Then the effective filling is 
u* = 1 {p = 1) OT u* = —1 {p = 2) and the Hartree-Fock contributions are the same in 
the two descriptions. The RPA collective modes spectra look very similar - Figures 3-4 
(the agreement is exact at g — ^ 0) but the values of interaction energy differ a lot (Table 
III and Table IV). 

In summary the agreement between the RPA interaction energies and the exact 
results are far from the expected one and the analysis needs an extension by including 
three-body contributions (three-body density- density correlation function |T^) which 
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seems to be very complicated. An alternative approach may be obtained within the new 
formalism developed by Shankar and Murthy (which has a direct relation with the 
Laughlin trial wave function). In Ref. |2T| they showed how to calculate the energy gaps 



and their results agree reasonably with numerical results |2 
4. Conclusions 

The values of Coulomb interaction energies for the 2D electron system in the region 
of FQHE are calculated within the Chern- Simons theory in the RPA for several frac- 
tional fillings. The results are obtained in the limit — ^ (i. e. , when 
Coulomb interaction is very small comparing to the separation between Landau levels) 
and compared with the exact diagonalization results (the results for few particle systems 
extrapolated for infinite systems). The best agreement is found for fractions going down 
from 1/2 (3/7, 2/5, 1/3), for the best 1/3 result the difference between the RPA and the 
exact results is of order of 20% (of the exact value). A qualitative difference is obtained 
for u = 1 p = 1 CF description when a finite RPA correlation energy is found. Also 
the result for i^ = l/3p = 2is very different from = 1/3 result obtained for p = 1. 
Our analysis needs an extension to a higher order approximation including three-body 
contributions. An alternative approach may be obtained within the new formalism de- 
veloped by Shankar and Murthy [Q, their results for energy gaps |21| are in reasonable 
agreement with numerical results ^2 . 
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Figure 1: 

Collective modes for i/ = 1, = 1. 
Figure 2: 

Collective modes for the filling u — 1 given within the p — 1 CF description 
r. = l. 

Figure 3: 

Collective modes ior u — 1/3, Vg — 1. 
Figure 4: 

Collective modes for the filling v = 1/3 given within the p = 2 CF description 
r. = l. 

Figure 5: 

Collective modes for i/ = 3/7, = 1. 
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H-F 


exact 


1 


-0.627 


-0.627 


2/3 


-0.497 


-0.519 


5/11 


-0.406 


-0.451* 


4/9 


-0.402 


-0.447* 


3/7 


-0.396 


-0.443 


2/5 


-0.385 


-0.433 


1/3 


-0.362 


-0.412 


1/5 


-0.280 


-0.328 



Table I: The Hartree-Fock and exact interaction energies (per particle) in respective units 
of The "exact" results are taken from Refs. Q where the results (in spherical 
systems) for few particles (A^ < 12) were extrapolated for infinite systems (A^ oo). 
Two numbers with stars are obtained within the Jain CF approach 0. The "exact" 2/3 
result is found via particle-hole symmetry 
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V 






W 








= 10-4 


Ts = 10-5 








h — in 


h — 1 ^ 
Av — lO 


h — on 


h — in 


Av — lO 


h — on 

Av — ZU 


3/7 


3 


1 


— U.OOl 


— u.ozo 


— u.ozu 


— U.OOl 


— u.ozo 


n f^on 
— u.ozu 


2/5 


2 


1 


-0.360 


-0.357 


-0.355 


-0.360 


-0.357 


-0.355 


1/3 


1 


1 


-0.247 


-0.247 


-0.246 


-0.247 


-0.247 


-0.246 


1/5 


1 


2 


-0.548 


-0.544 


-0.542 


-0.548 


-0.544 


-0.542 


1 


-1 


1 


-0.230 


-0.230 


-0.230 


-0.230 


-0.230 


-0.230 


1/3 


-1 


2 


-0.585 


-0.581 


-0.580 


-0.585 


-0.581 


-0.580 


2/3 


-2 


1 


-0.358 


-0.354 


-0.352 


-0.358 


-0.354 


-0.352 



Table II. The RPA correlation energies in units of — 177. For the A; = 10 case of z/ = 3/7 
the summation goes over 15 modes. For negative effective fillings we used instead u* > 
p < (1/ < 0). 



u 


u* 


P 


N 


N 


N ' N 


exact diagonalization 


1 







-0.627 





-0.627 


-0.627 


3/7 


3 


1 


-0.396 


-0.197 


-0.593 


-0.443 


2/5 


2 


1 


-0.385 


-0.159 


-0.544 


-0.433 


1/3 


1 


1 


-0.362 


-0.142 


-0.504 


-0.412 


1/5 


1 


2 


-0.280 


-0.242 


-0.522 


-0.328 



Table III. The interaction energies in respective units of (note the change of units 
with respect to Table II). The RPA values in the fifth column are taken from A; = 20 
results of Table II. 
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V 


V* 




i V 


/\7 

i V 


AT ' AT 


exact diagonalization 


1 


-1 


1 


-0.627 


-0.230 


-0.857 


-0.627 


1/3 


-1 


2 


-0.362 


-0.335 


-0.697 


-0.412 


2/3 


-2 


1 


-0.497 


-0.203 


-0.700 


-0.519 



Table IV. The interaction energies in units of obtained for negative effective fillings. 
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Figure 1 




Figure 2 




Figure 3 




Figure 4 




Figure 5 



